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Abstract 

We perform a fully quantum mechanical numerical calculation for the problem 
of a single electron (or excitation) propagating in a iV-site one- dimensional chain 
in the presence of a single Holstein impurity. We compute the long-time averaged 
probability for finding the electron on the impurity site as a function of the nonlin- 
earity parameter, defined in terms of the electron-phonon coupling strength and the 
oscillator frequency. The results, in the intermediate nonlinearity parameter range, 
differ substantially from the ones obtained through the use of the discrete nonlinear 
Schrodinger equation, even in the high-frequency regime. 
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Interest concerning issues related to the interplay between nonlinearity and dis- 
order in discrete systems lead us to consider some time ago, the problem of a single 
nonlinear (cubic power) impurity embedded in a discrete, linear, one-dimensional lat- 
tice. By using the lattice Green functions in a selfconsistent manner we were able 
to derive analytically all the bound state properties and the transmission coefficient 
through the impurity]!], |2|. We found an impurity nonlinearity threshold below which, 
there is no bound state, in marked contrast with the linear impurity case||. The dy- 
namics of the problem showed the presence of a selftrapping transition, when the 
impurity nonlinearity exceeded a threshold value H f§, ||. We also considered the case 
of a generalized impurity and, by using either an appropriate ansatefd or the Green 
function formalism were able to construct analytically a phase diagram in param- 
eter space showing regions with different number of bound states. Extensions of the 
nonlinear impurity problem to a dimerized chain and to higher dimensions have been 
recently investigated ||. 

In all of the above studies, the evolution of a quasiparticle propagating in the 
lattice is given by the discrete nonlinear Schrodinger equation (DNLS): 

dC 

i-^ = nCn-i + C n+ i) - x|C„r _1 C n 5 n0 (1) 
dt 

where C n (t) denotes the probability amplitude for the quasiparticle to be at site n, 
V is the nearest-neighbor transfer energy integral, \ is the nonlinearity parameter 
and a is the nonlinearity exponent. In the Holstein impurity context, the DNLS (for 
the cubic nonlinearity, a = 3) describes the time evolution of an electron which is 
strongly coupled to the vibrational degrees of freedom of a molecule at the impurity 
site. In this case \ = g 2 /uj where g is the electron-phonon coupling and u is the 
phonon frequency. 

However, several doubts have been raised concerning the general applicability of 
the semiclassical approximation on which the DNLS is based to the fully quantum 
mechanical Holstein problem[9 - 13]. This has motivated us to calculate the effect 
of quantum fluctuations in these impurity systems. In the present letter we perform 
a fully quantum mechanical numerical calculation for a single vibrational impurity 
coupled to an electron (or excitation) at a given site (j = 0) of a finite chain with 
periodic boundary conditions and compare the results to the DNLS predictions. The 
Hamiltonian is 

JV-l 

H = V Y, {CjCj+i + C] +1 Cj} + uata + g{a + a f ) CjO, (2) 

3=0 

where Cj(Cj) creates (destroys) an electron on the jth-site, aJ(a) is the usual phonon 
creation (destruction) operator, uj is the (Einstein) oscillator frequency and g is the 
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electron-phonon coupling at the impurity site. The most general eigenfunction has 
the form 

N-l oo 

m = EE«i"V> ® b> (3) 

where \v) = {a)) v / y/v\ \fy P h and \j) = Cj |0) e . The states \0) p h and |0) e are the 
phonon and electron vacuum, respectively. After inserting Eq.fl3|) into the eigenvalue 
equation H\^e) — E \^?e), we obtain 

N-l oo N-l oo 

ED^-^l 4 )c » = t/ EEl 4-± + ]<?» 

j=0 u=0 j=0 u=0 

oo 

+ gel £ [ «!f +1) vVTI + at 1) V^]W) (4) 

i/=0 

The use of the transfer matrix formalism and periodic boundary conditions allows us 
to express all the aq' in terms of the viz., 

4 = I sin C^) + sin ((^ - 3)hu) ] (5) 

where the wavevector k v is defined by A u = (E—u v)/2V = cos(fc„) and the eigenvalue 
equation reduces to 

(#/V) sin(iVA; l/ ) {c^ +1) VlTTl + ajp 1 ^} = 2«M sin(^) [ cos(Nk u ) - 1 ] (6) 
When > 1, k v — > i K u + {0,7r} where the left (right) choice corresponds 



to A v > 1 (A, < -1) and where K v is given by K v = log(\A u \ + JAl - 1). The 
trigonometrical equations fl5|), (§) change to a hyperbolic form: 



«f = ^jfjfjQ [ (±1)" sinh(j^) + sinh((iV- j)^) ] (±1) J (7) 



sinh(Ar^){a^ +1) ^^7+l + a { o~ 1] V^} = ±2 sinh(if ly ){cosh(^ i/ )-(±l) iV } 

(8) 

where the ± sign corresponds to the sign of A v . Equations @ and (§) constitute an 
infinite set of "three-terms" recursion relations of the form: 

C a {u ~ 1] + b u a {v) + a u a {u+1) = (9) 

The eigenvalue equation associated to this (infinite) tridiagonal equation can be solved 
by a standard method. However, some care must be taken in order to avoid numerical 
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instabilities: we iterate Eq.fl^) 
"down" from u mn7 . down to L, 



from two directions, "up" from v = up to L and 
where v max is an appropriate phonon number cutoff. 
The eigenvalue equation is then obtained by imposing that both iterations yield the 
same value for a^ +1 ^ /a^ at the "gluing" point v — L, chosen to ensure the numerical 
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stability of the iteration process, 
monitored through several tests: 
(a) Evaluation of 



The precision of the numerical computation is 
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(L+l) [ ) (L+l) 
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(10) 



(b) Orthogonality relations: for instance, J2u a o^! = where the index fj, runs 
over all the eigenstates. 

(c) A general test based on the Hamiltonian: If A is an arbitrary operator, then the 
mean value of Q = [A, H] over any eigenstate of H is identically zero. For instance, 
if we take A = a, the phonon destruction operator, the identity = (Q) leads to the 
condition u < a > + g < CqCq >= 0. 

Once in possession of all the eigenenergies and the amplitudes oq^ we are 
ready to compute dynamical observables. As the initial condition for the electron 
we use the one that places it completely on the impurity site at t — and focus on 
the probability Po(t) for finding it there at an arbitrary time t later. For the phonon 
part we use two different initial conditions: an undisplaced oscillator (zero phonons 
present) and a naturally "relaxed" oscillator (a coherent state). Due to absence of 

degeneration, the long-time average probability < P >= liniT-»oo 0-/T) Jq Po(t) dt 

0) 



can then be expressed in terms of the a ll as 



Initial state: |^(0) >= \j = > <g> ^p n \n > 



n=Q 



and result 



where 



Pn 



i f-g 



n\ \ to 



<Po>= EE 



a 



0, fi 



(n) t2 



n=0 



undisplaced oscillator 



exp[— (g/id) 2 /2] coherent state case 



(11) 



(12) 



We have examined the behavior of Eq.(fTT|) as a function of the nonlinearity pa- 
rameter \jV = g 2 /ujV, for different oscillator frequencies (u oc 1/^mass), for chains 
with = 2, 3, 4 and 5, the chain corresponding to iV = 2 being the only nonperiodic 
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one. We show the results for N = 2 and N = 4 in figures 1-4 where we also compare 
with the quantity < P > evaluated though the DNLS equation. We observe signifi- 
cant departure from the DNLS prediction for intermediate nonlinearity, even for very 
large u. In the case of an (initially) undisplaced oscillator, while < Pq > increases 
with x on average, it also shows oscillatory behavior with \ for large oscillator mass. 
For smaller masses, the oscillatory behavior decreases setting into a monotonic curve 
for the smallest mass examined. This case, which is closer to the antiadiabatic limit 
where the DNLS equation should work, still differs significantly from the DNLS pre- 
diction. For the case of an (initially) displaced oscillator, < Pq > increases always 
monotonically with the nonlinearity parameter. In both the oscillator mass 

decreases, the curves converge to a linear-like impurity curve Jl|, [|, remaining signif- 
icantly away from the DNLS prediction, down to the lowest mass considered here. 
For io = 100V the curve is almost undistinguishable from the one corresponding to a 
static, linear impurity of strength x- We note, in particular, the dramatic difference 
in the slopes of the quantum mechanical versus the DNLS curve in the small coupling 
limit. The quantum mechanical impurity results approach the DNLS impurity results 
only in the very large coupling limit \ 00 • Similar results are found for the N = 3 
and N = 5 case. 

The equivalence of model (|2|) with a static impurity can be understood as follows: 
Let us perform the canonical transformation H —± H = U H U^ 1 with U = exp[ A(a — 
at) Co Co]. With the choice A = —g/uj, the transformed Hamiltonian becomes: 

H = Y,[ VjCjCj+i + h.c. }+uj<x ] a- (g 2 /u) C f C (13) 



j 



where 



V exp[ (g/w)(a-at) } j = 
V 5 ={ V exp[ -{g/u){a - at) ] j = -1 (14) 

V jV0,-i 

In the limit u — ► oo with fixed x = 9 2 / u 

H -»■ J2 [ V C}C j+1 + h.c. ]+u a f a - {g 2 /u) C f C (15) 
j 

which is the Hamiltonian for a chain with a single static, linear impurity. Therefore, 
in the limit uj —>■ oo (u ~ 100 V is good enough) with fixed x — 9 2 l^i our Hamiltonian 
becomes equivalent to the one of a single, static, linear impurity. The DNLS limit is 
never achieved. 

We have directly shown in the present work that the problem of a single vibrational 
impurity embedded in a finite linear chain, when treated in a completely quantum- 
mechanical framework, does not approach the predictions based on the DNLS. We 
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conclude that, for the Holstein impurity problem the DNLS fails to describe cor- 
rectly the selftrapping process. Our results agree with recent studies[[L3|] on a coupled 
quasiparticle-boson system which also finds that the DNLS has a limited range of 
validity in the context of this problem. It would be interesting to test the DNLS in 
a system possessing translational symmetry where the enhancement of trapping due 
to the single impurity character of the present model would be absent. Work in that 
direction is currently under way by our group. 

Finally, It is worthwhile pointing out that the general solutions for DNLS-like 
impurity problems that have been obtained earlier retain their validity and usefulness 
in classical problems where DNLS is valid. Examples of these problems are typically 
furnished from nonlinear optical applications. 

Two of the authors (M. I. M. and J. A. R.) acknowledge partial support from 
FONDECYT grant 1950655 while the third one (G. P. T.) acknowledges support from 
UENEA grant 95-i?A-115 of the General Secretariat of Research and Technology of 
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FIGURE CAPTIONS 



Figure 1: N = 2 chain: Long-time average probability of finding the electron (exci- 
tation) on the impurity versus the nonlinearity parameter xjV = g 2 /uV, for different 
oscillator frequencies (masses). The oscillator is initially undisplaced. 

Figure 2: Same as in figure 1, but for the case of an initially displaced oscillator. 

Figure 3: N = A chain: Long-time average probability of finding the electron (exci- 
tation) on the impurity versus the nonlinearity parameter xjV = g 2 /uV, for different 
oscillator frequencies (masses). The oscillator is initially undisplaced. 

Figure 4: Same as in figure 3, but for the case of an initially displaced oscillator. 
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